Transport properties of two-dimensional electron systems on silicon (111) surfaces 
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We theoretically study transport properties of a two-dimensional electron system on a hydrogen- 
passivated Si(lll) surface in the field-effect-transistor (FET) configuration. We calculate the density 
and temperature dependent mobility and resistivity for the recently fabricated Si(lll)-vacuum FET 
by using a semiclassical Boltzmann theory including screened charged impurity scattering. We find 
reasonable agreement with the corresponding experimental transport properties, indicating that the 
screened disorder potential from random charged impurities is the main scattering mechanism. We 
also find that the theoretical results with the valley degeneracy g v = 2 give much better agreement 
with experiment than the g v = 6 situation indicating that the usual bulk six-valley degeneracy of 
Si is lifted in this system. 
PACS Number : 73.40. Qv, 72.10.-d, 73.40.-c 



Recently a high mobility two dimensional electron sys- 
tem (2DES) on Si(lll) surface has been fabricated [1,2]. 
The fabricated Si- vacuum field-effect transistor (FET), 
in which an electric field is applied through an encapsu- 
lated vacuum cavity instead of the usual SiC>2 as in a Si 
MOSFET, allows the 2DES to be gated on a hydrogen- 
passivated Si(lll) surface. The strong scattering poten- 
tial associated with the disordered Si-SiC>2 interface in 
the usual Si MOSFET geometry [3] being absent at the 
H-passivated Si(lll)-vacuum interface, the carrier mobil- 
ity is expected to be higher in this new 2D system com- 
pared with the extensively-studied [3] Si(lll)-Si02 2D 
MOSFET system. In the usual Si-Si0 2 MOSFETs the 
carrier mobility is limited by charged impurity scattering 
at low carrier densities, but at higher carrier densities, 
when the 2DES is pushed very close to the Si-Si02 inter- 
face by the steep self-consistent confining potential aris- 
ing from the carriers themselves, surface roughness and 
defect scattering associated with the Si-Si02 interface be- 
comes dominant, reducing carrier mobility with increas- 
ing density. The 2D mobility in Si-Si0 2 MOSFETs there- 
fore has a peak value at some intermediate carrier density 
(~ 10 12 cm -2 ) with the mobility limited by charged im- 
purity scattering (interface roughness scattering) below 
(above) the characteristic density where the 2D mobil- 
ity attains its peak value. Potential advantages of the 
newly created [1,2] Si- vacuum FET structure are that the 
absence of interface roughness/dangling bond scattering 
leads to (1) very high 2D mobilities not achievable in Si- 
Si0 2 MOSFET structures; and (2) the mobility should, 
in principle, be a monotonically increasing function of 
2D carrier density, as, for example, in 2D undoped GaAs 
HIGFET structures [4]. In the Si- vacuum FET fabri- 
cated by Eng et al. [1,2] the H-Si(lll) surface reduces 
surface scattering, thereby increasing electron mobility. 
The measured [1] peak electron mobility on H-Si(lll)- 
vacuum FET reaches ~ 8000 cm 2 /Vs at 4.2K, which 
is the highest electron mobility recorded on a Si(lll) 
surface compared to the previous peak electron mobility 



measurements in Si(lll) MOSFETs (~ 2500 cm 2 /Vs) 
[5-7]. With this enhanced electron mobility, the Si- 
vacuum FET has potential development of atomic-scale 
devices, where electrons are coupled to specific molecules 
or single atoms positioned on the H-passivated Si surface. 
Indeed, with further improvement in the surface passiva- 
tion techniques developed in refs. [1,2], the experimental 
mobilities should go up much higher than the currently 
published [1] record of ~ 8000 cm 2 /Vs. 

In contrast to the well-studied electron inversion lay- 
ers in Si(100) MOSFETs [3], the detailed investigation 
of electronic transport properties of 2DES on Si(lll) 
surface has not been reported in the literature. Since 
the high mobility FET on Si(lll) surface is now avail- 
able, a detailed study of electronic transport properties 
is necessary. The purpose of this paper is to calculate 
the density and temperature dependent electron mobil- 
ity of the 2DES on Si(lll) surface. We also calculate 
the 2D temperature-dependent resistivity in connection 
with the subject of 2D metal-insulator transition and 
the 2D metallic phase which has been studied for most 
other high mobility 2D systems [8]. The measured tem- 
perature dependent resistivity in Si-vacuum FET shows 
strong metallic behavior [2] (i.e. the resistivity increases 
with temperature) with no sign of insulating behavior 
down to a density n = 3.7 x 10 n cm~ 2 (which is the low- 
est achieved experimental density so far). The density 
dependence of mobility in the metallic regime of the Si- 
vacuum FET may give valuable information about the 
scattering mechanisms operational in limiting 2D carrier 
transport, since different scattering sources give rise to 
different density dependence in the resistivity. Thus, 
a detailed study of the density and temperature de- 
pendence of 2D resistivity in high-mobility 2D Si(lll)- 
vacuum FET could be useful not only for a better fun- 
damental understanding of the limiting transport mecha- 
nisms in this potentially important new system, but also 
for providing ideas about how to further increase the 2D 
mobility leading to better devices. 
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In order to calculate the low temperature 2D electron 
mobility we have carried out a microscopic transport 
calculation using the Boltzmann transport theory [8]. 
We calculate the mobility in the presence of randomly 
distributed Coulomb impurity charges near the surface 
and the bulk-acceptors in the depletion layer with the 
electron-impurity interaction being screened by the 2D 
electron gas in the random phase approximation (RPA). 
The screened Coulomb scattering is the only important 
scattering mechanism in our calculation. There are ad- 
ditional interface-scattering mechanisms unrelated to the 
Coulomb centers (e.g. surface roughness scattering), but 
such interface scattering would become more pronounced 
at higher densities, as the inversion layer is drawn closer 
to the interface. In general, the surface roughness scatter- 
ing gives rise to the decrease of the mobility at high den- 
sities. The observed mobility [1] shows that the mobility 
increases monotonically with electron density indicating 
the apparent absence of interface scattering, as expected 
for a Si- vacuum interface. In fact, as argued above in this 
paper, the surface roughness scattering should be negligi- 
ble in the Si- vacuum FET systems, allowing us to ignore 
this process in our theory. We also neglect all phonon 
scattering effects mainly because our theoretical estimate 
shows phonon scattering to be negligible for Si-vacuum 
FET structures in the T < 1QK regime of interest to 
us. We note that our theoretical approach has provided 
a reasonable description of the transport properties for 
several 2D systems [8]. 

In the Boltzmann theory the mobility is given by 
ji = e(T)/m with m as the carrier effective mass, and (r) 
the energy averaged finite temperature scattering time: 
(r) = jde k e k r{e k ) (~ / f de k e k wh cre f(e k ) 

is the Fermi distribution function, f(e k ) = {1 + exp[(e£ — 
/j)]/^T} _1 with n(T, n) as the finite temperature chem- 
ical potential determined self-consistently. The energy 
dependent scattering time r(e k ) for our model of ran- 
domly distributed impurity charge centers is given in the 
leading-order theory by 

1 2tt ^ f (o) f d 2 k' v{q,z) 2 

x(l-cos0)$(e k -e k /), (1) 

where nf" 1 is the concentration of the a-th kind of charged 
impurity center, q = |k — k'|, 9 = # kk ' is the scattering 
angle between the scattering in- and out- wave vectors 
k and k', e k = h 2 k 2 /2m, v(q,z) is the 2D Coulomb 
interaction between an electron and an impurity at z, 
and e(q) = s(q; fi, T) is the 2D finite temperature static 
RPA dielectric (screening) function [3]. In calculating 
the Coulomb interaction and the RPA dielectric function 
in Eq. (1) we take into account subband quantization 
effects in the inversion layer through the lowest subband 
variational wavefunction [3]. From calculated mobility 



we get the resistivity given by p = l/ne/j,. 

Throughout this paper we use the following material 
parameters: valley degeneracy g v — 2 (or g v — 6 as the 
case may be), an effective mass m — 0.35m e (correspond- 
ing to the Si(lll) surface) where m e is the free electron 
mass, usi = 11.7 and K vac = 1, the dielectric constant 
of Si and vacuum, respectively. Even though the valley 
degeneracy (g v ) on the Si(100) surface is well-known to 
be g v = 2, the valley degeneracy on the Si(lll) surface 
has been a puzzle for a long time [3]. Since Si has six 
conduction-band minima located along the (100) direc- 
tions, the (111) surface should have six equivalent valleys 
(i.e. g v — 6) within the simple effective mass approxi- 
mation. However, magnctotransport measurements have 
found the valley degeneracy to be both 2 and 6 [5-7,3], 
with most experiments finding the lower valley degener- 
acy of 2. In our calculation we get much better agreement 
with the experimental data of refs. [1,2] when we use the 
valley degeneracy g v = 2 instead of g v = 6, indicating 
that the valley degeneracy is reduced from 6 to 2 in the 
Si(lll)-vacuum FET structures of rcf. [1]. 

In Fig. 1 we show the calculated mobility as a function 
of electron density for different temperatures. The mo- 
bilities for Si(lll)-vacuum FET are shown in Fig. 1(a) 
with valley degeneracy g v = 2 and Fig. 1(b) with g v = 6. 
In Fig. 1(c) we show the mobility of a Si(lll) MOSFET 
system with g v = 2. For the Si-Si0 2 MOSFET calcu- 
lations we use Ksio 2 — 3.9 with the other parameters 
being the same as in the Si-vacuum system. In Fig. 1 
the mobilities are calculated with fixed impurity density 
rii = 7 x 10 11 cm~ 2 (the impurity centers are distributed 
completely at random in a 2D plane at d — 50A) and 
the ionized impurity density ti^d = 4 x 10 16 cm~ 3 in 
the depletion layer of the p-type Si. Here d indicates the 
distance between the impurity charge centers and the in- 
version layer. This is intended to represent the average 
distance between an electron in the inversion layer and 
the impurity near the Si-vacuum interface. The calcu- 
lated mobilities depend on the density approximately as 
/i ~ n a for the calculated density range. We emphasize 
that this is not a strict power law since the exponent a 
depends weakly on the density. Over the experimental 
density range, we find that for g v = 2 [Fig. 1(a)] the 
exponent a in /i oc n a is about 0.7 at T=0K and in- 
creases as temperature increases. At T = AK we have 
a w 1.2 which is consistent with experiment [1,2]. For 
g v = 6 [Fig. 1(b)] we have a w 0.5 - 1.5. Note that the 
Si(lll) MOSFET system [Fig. 1(c)] has a « 0.3 - 0.6 in 
the given temperature range. [The characteristic separa- 
tion "d" between the charged impurities and the 2DES 
as well as the random impurity densities rij and n^o are 
of course not precisely known for the experimental sam- 
ples, but the parameters we use are quite typical, and 
give decent agreement with the experimental results of 
ref. [1] for g v — 2, i.e. Fig. 1(a).] Our mobility calcu- 
lation indicates that the screened potential from random 
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FIG. 1. Calculated electron mobility in the Si(lll)-vacuum 
FET system for various temperatures, T = 0, 1, 2, 3, 4K 
(top to bottom) with valley degeneracy (a) g v = 2 and (b) 
g v — 6. Fig. 1(c) shows the calculated electron mobility 
for a Si(lll)-Si0 2 MOSFET with g v = 2. The calculated 
mobilities depend on the density approximately as /i ~ n a 
with a « 0.7 at T = and a « 1.2 at T = 4i(T for g„ = 2 
[Fig. 1(a)], and a w 0.5 - 1.5 for g v = 6 [Fig. 1(b)]. For 
Si(lll) MOSFET [Fig. 1(c)] we have a ^ 0.3 - 0.6. Note 
that surface roughness scattering, known to be important for 
Si-Si02 MOSFETs at higher carrier densities, is neglected in 
the calculation. 
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FIG. 2. Calculated resistivities for various hole densities 
n =3.7, 4.0, 4.5, 4.9, 5.36 x 10 11 cm -2 (from top to bot- 
tom) with the same parameters of Fig. 1. The dashed 
lines show the calculated resistivity for density n =4.5, 4.9, 
5.36x 10 n cm~ 2 with the valley degeneracy g v = 6. 

charged impurities is the main scattering source and the 
Si(lll)-vacuum FET has a valley degeneracy g v = 2. 
Since there are no adjustable parameters in our calcu- 
lation except the impurity parameters, the agreement is 
quite good. The calculated mobility obviously depends 
on the distance (d) of the charged impurity centers from 
the 2D layer. Generally the location of the charged im- 
purities in the FET is not known. Thus, the choice of 
d = 50 A is reasonable because our calculated mobil- 
ity has the same dependence on n as does the measured 
mobility [1,2]. 

In Fig. 2 we show the calculated temperature depen- 
dent resistivity using the same parameters as in Fig. 1. 
The calculated resistivity increases as the temperature 
increases (i.e. metallic behavior), and shows an approx- 
imately linear temperature dependence in the T < 5K 
regime as observed in the experiment [2]. The tem- 
perature dependence of p(T) becomes stronger with de- 
creasing density just as the experimental observation in 
Si(lll)-vacuum FET [2]. This strong linear temperature 
dependence of p(T) in this system arises [8,9] due to the 
peculiar nature of the 2D screening function which has 
a cusp at wave vector q = 2kp. Since 2fcp scattering is 
the most dominant resistive scattering at low tempera- 
tures, thermal smearing of this 2kp cusp leads to strong 
temperature dependence of resistivity. Thus, tempera- 
ture dependent screening of charged impurity scattering 
in 2D layers leads to a strongly temperature dependent 
effective disorder controlling p(T, n) at low temperatures 
and densities. (The screening effect decreases with in- 
creasing temperature giving rise to increasing effective 
disorder with increasing T, and hence increasing p(T) 
with T.) Note that the overall resistivity scale depends 
on the unknown impurity densities, but the qualitative 
trends in p(T, n) arise from basic aspects of the underly- 
ing scattering mechanisms. 
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We pointed out earlier [8] that the requirements for 
the observation of a large temperature-induced change 
in resistivity is the strong screening condition, qtf 3> 
2k p, where qxF and Uf are the Thomas- Fermi screening 
wave vector and Fermi wave vector, respectively. Since 
9tf/^f oc g^ 2 we expect the observed metallicity to be 
much stronger for g v = 6 than for g v = 2. The dashed 
lines in the Fig. 2 show the calculated resistivity for 
n =4.5, 4.9, 5.36xlO n cm~ 2 with g v = 6. We find, in the 
range of of — 5K, roughly a factor of three change in 
p{T) with g v — 6 whereas more like a 50% temperature 
induced change in p(T) with g v = 2, which is consis- 
tent with the experimental observation. Thus, the trans- 
port measurement in Si(lll)-vacuum FET supports the 
Si(lll) surface having a valley degeneracy 2. 

Finally, we briefly discuss the possible origin of the val- 
ley degeneracy puzzle in this system, i.e. why g v = 2 in 
the Si(lll)-vacuum FET system rather than the usual ef- 
fective mass prediction of g v = 6. This is, in fact, an old 
problem [3] much discussed in the context of Si(lll)-Si02 
MOSFETs. It is generally accepted [5,6] that the random 
strain associated with the Si surface lifts the valley de- 
generacy lowing two of the six valleys in energy. Only 
specially processed Si(lll)-Si02 samples with rather low 
mobilities were ever found to have g v = 6 [6] with most 
Si(lll)-Si0 2 MOSFETs having g v '= 2. We believe the 
same kind of one electron physics to be operational in 
lifting the valley degeneracy from g v — 6 to g v — 2 in the 
Si(lll)-vacuum system as well although understanding 
the details of this valley degeneracy puzzle is beyond the 
scope of our current work where we are interested in the 
transport properties. We do, however, mention that this 
valley degeneracy lifting is unlikely to be a many-body 
exchange instability, which would occur even as a mat- 
ter of principle at a much lower critical density (<~ 10 11 
cm~ 2 ) with a critical temperature of only IK or so. At 
the experimental density and temperature range of ref. 
[1] , a many-body effect induced valley degeneracy lifting 
is extremely unlikely [10,11] 

In conclusion, we calculate the the density and tem- 
perature dependent mobility and resistivity for the re- 
cently fabricated Si(lll)-vacuum field effect transistor 
by using a semiclassical Boltzmann theory including 
screened charged impurity scattering. We find reasonable 
agreement with the corresponding experimental trans- 
port properties for the (111) system indicating that the 
screened potential from random charged impurities is the 
main scattering source. In comparing our calculation 
with the existing experimental data [1,2] we find that 
the calculated results with the valley degeneracy g v = 2 
give better agreement with experiments than g v = 6. 
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